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ABSTRACT 


As  an  aid  Lo  better  understanding  ELF  propagation  in  the  earth- 
ionosphere  waveguide,  full  wave  ELF  modal  height  gains  are  analyzed  as 
regards  regions  of  reflection  and  absorption.  Frequencies  in  the 
Schumann  and  Seafarer  bands  are  discussed  for  representative  day/night, 
ambient  and  disturbed  ionospheres.  Enhancement  of  the  earth  ionosphere 
waveguide  attenuation  due  to  sporadic  E layering  is  discussed  and  the 
origin  of  the  enhancement  clarified. 


1.  INTRODUCTION 

In  the  field  of  extremely  low  frequency  and  very  low  frequency  (ELF), 
(VLF)  propagation  in  the  earth-ionosphere  waveguide,  much  attention  has 
been  given  in  recent  years  to  the  problem  of  solving  the  modal  equation 
for  ionospheres  having  a gradual  continuous  increase  in  electron  and  ion 
densities  with  height  at  the  lower  boundary  while  simultaneously  making 
allowance  for  anisotropy  resulting  from  the  geomagnetic  field.  Most  of 
these  solutions  are  obtained  numerically  and  various  approaches  are 
documented  in  the  works  of  Budden^’^,  Wait^,  Galejs^  and  Krasnushkin^.  A 
by-product  of  numerical  solutions  is  much  quantitative  information  con- 
cerning critical  ionospheric  height  regions  as  regards  reflection  and 
absorption  processes.  Particularly  as  regards  ambient  ionospheres,  where 
the  geomagnetic  field  plays  a role,  this  information  has  been  to  a great 
extent  ignored  because  of  the  intense  iriterest  in  the  normal  mode  proper- 
ties themselves  (i.e.  phase  velocities  and  attenuations)  as  well  as 
their  excitatiorrs.  This  is  unfortunate  since  detailed  full  wave  analysis 
of  nomal  mode  reflection  and  absorption  can  provide  valuable  Insights 
into  puzzling  features  of  the  modal  properties.  As  an  example,  from 
full  wave  outputs  for  ELF  nomal  mode  attenuation  rates  it  has  been 
observed  that  the  presence  of  a sporadic  E layer  at  nighttime  could  in- 
crease attenuation  from  an  order  of  magnitude  of  1 dB/Mm  to  about  10  dB/Mm 
in  the  Seafarer  band^»^.  But  without  detailed  full  wave  analysis  of 
nomal  mods  reflection  and  absorption,  the  underlying  mechanisms  must 
remain  somewhat  obscure.  Mor'eover,  the  attenuation  enhancement  cited 
above  may  play  a significant  v'ole  in  the  relatively  short  path  (1600  km) 
6-0  do  fades  reported  In  connection  with  emissions  from  the  Wisconsin  Test 
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Facility  ’ . Finally,  a detailed  full  wave  analysis  of  normal  mode  reflec- 
tion and  absorption  can  provide  coi.'^tructive  guidance  in  the  problem  area 
of  profile  fitting  to  radio  frequency  data  in  the  ELF  band.  Thus,  to 
illustrate  some  of  the  information  available  from  full  wave  programs 
as  well  as  to  provide  a better  understanding  of  ELF  propagation  in  the 
earth-ionosphere  waveguide,  full  wave  ELF  modal  height  gains  are  analyzed 
in  this  report  for  a variety  of  ionospheres  as  regards  regions  of  reflec- 
tion and  absorption.  Results  will  be  presented  for  frequencies  between 
the  Schumann  band  and  about  100  Hz  for  representative  day/night  ambient 
and  disturbed  ionospheres. 

In  the  following  section  notation  for  some  of  the  more  important 
symbols  is  given.  In  section  3 relevant  formulas  utilizing  the  full  wave 
fields  will  be  summarized,  in  section  4 results  are  presented  and  in 
section  5 conclusions  are  summarized.  Appendices  A and  B contain  deriva- 
tions of  formulas  useful  for  understanding  features  of  the  full  wave 
absorption  results  and  Appendix  C contains  the  derivation  for  a simplified 
sporadic  E layer  model  which  gives  insight  into  the  origin  of  attenuation 
rate  resonances  associated  with  such  layering. 

The  full  wave  fields  program  utilized  in  this  study,  apart  from 
modifications  effected  to  extract  field  informiation  in  the  most  compre- 
hensible form  for  our  purposes,  was  developed  by  Mr.  C.  H.  Shellman  and 
is  a Runge-Kutta  variant  of  the  original  fields  program  developed  by 
Mr.  R.  Smith^O.  The  latter  was  based  on  the  Inoue-Horowitz^^  and  Pitte- 
way^2  papers.  The  earth-ionosphere  waveguide  eigenangles  which  serve  as 
input  to  the  full  wave  fields  program  and  modal  attenuation  rates  pre- 
sented In  section  4 are  based  on  Budden's  formal ism^*^  as  discussed  by 
Pappert  and  Moler^^. 
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2.  NOTATION 


Unless  otherwise  specified  MKS  units  are  assumed  throughout  the  study 
reported  here  and  some  of  the  more  important  symbols  are  the  following: 


X,  y,  z 


E 

H 

— k — i 

3C  = ZqH 
P 

Tt 

Iq 

Ei 

Hi 

3q 

Pi 

»i 

* 

Iiu 

Re 

Kq 

(i> 


u 

t 


Cartesian  coordinates;  z is  measured  vertically  up- 
wards into  the  ionosphere  and  x is  measured  horizon- 
tally with  x-z  the  plane  of  incidence. 

Radio  wave  electric  field  intensity. 

Radio  wave  magnetic  field  intensity. 

Normalized  radio  wave  magnetic  field  intensity. 
Polarization  vector. 

Poynting  flux  vector. 

Characteristic  impedance  of  free  space, 
i = 1,  2,  3;  X,  y,  z components  respectively  of  L 

i » 1,  2,  3;  X,  y,  z components  respectively  of  Tf. 

i = 1,  2*  3;  X*  y,  z components  respectively  of  K. 

i » 1,  2,  3;  X,  y,  z components  respectively  of*?, 

i « I,  2,  3i  X,  y,  z co<iiponents  respectively  of  n. 

Complex  conjugate. 

Imaginary  part  of. 

Real  part  of. 

Free  space  permittivity. 

2ifx  frequency  of  rf  wave. 

i » I,  2,  3;  j =*  I,  2,  3;  xx,  xy,  xz,  etc. . components 
of  susceptibility  tensor. 

Relative  heating  loss  per  unit  volume. 

Time. 
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2Tr/Ao  free  space  wavenumber  of  radio  wave. 

Free  space  wavelength  of  radio  wave. 

Plane  wave  angle  of  incidence  measured  from  z (can  be 

eigenangle) 

sin(e) 

cos(0) 

Collision  frequency  for  species  j. 

2t(X  cylotron  frequency  for  species  j (taken  to  be 
positive  for  all  species). 

2Ttx  plasma  frequency  for  species  j. 

Direction  cosine  between  geomagnetic  field  and  z axis. 
Speed  of  light  in  vacuum. 


^ ' 


3.  SUMMARY  OF  EQUATIONS 

Principal  outputs  of  full  wave  ELF  and  VLF  programs  such  as  documen- 
ted  by  Smith  are  the  electric  (E)  and  magnetic  (H)  fields  both  within  and 
without  the  ionosphere.  These  fields  could  be  identified,  as  is  often 
the  case,  with  vertically  or  horizontally  plane  polarized  waves  incident 
on  the  ionosphere  or  they  can  be  identified  with  generally  complex 
eigenangles  satisfying  the  fundamental  determinantal  mode  equation: 

F(0)=  11- (0)1=0  (1) 

with  I the  unit  matrix,  Rq  the  two  by  two  plane  wave  reflection  matrix 
associated  with  reflection  from  everything  above  the  level  D and  Rq 
the  plane  wave  reflection  matrix  associated  with  reflection  from  every- 
thing below  the  level  D.  Results  presented  in  the  next  section  are  in 
fact  for  fields  associated  with  the  zeros  of  equation  1 and  as  such  rep- 
resent normal  mode  reflection  and  absorption  properties.  Normaliza- 
tion of  the  modal  height  gain  functions  is  such  that  the  magnetic 
field  component  at  the  ground  is  1 volt/m.  In  the  remainder 
of  this  section  are  summarized  those  formulas  which  involve  the  full 
wave  fields  and  which  we  believe  are  most  useful  as  a diagnostic  tool  for 
understanding  fundamental  ELF  absorption  and  reflection  processes. 

a.  Absorption  and  attenuation. 

The  time  average  rate  of  conversion  of  rf  energy  into  heat  per  unit 


volume  is  given  by: 

j = p) 


Walts 


where  the  components  of  the  polarization  vector,  P,  are  given  by 


j=l 


»“  1,2,3 


The  time  average  of  the  Po.ynting  vector  is  given  by: 


1 „ “‘f  Watts 

= — Re  (E  X H 


and  is  generally  interpreted  as  the  time  averaged  energy  flux. 

A convenient  measure  of  the  integrated  effect  of  the  Joule 
heating  up  to  altitude  z is  in  terms  of  attenuation  units  given  by 


dB/1000  km 


where  h is  selected  high  enough  to  encapture  essentially  all  of 
the  energy  flux  in  the  x direction.  Although  lower  altitudes 
suffice*  results  of  Eq  (5)  given  in  the  following  section  are  generally 
based  on  an  h at  F region  altitudes  during  the  night  and  on  an  h of 
150  km  during  the  day. 


Results  showing  3a  as  a function  of  z will  be  given  in  the  following 
section.  Such  a display  will  show  the  buildup  of  the  absorptive  contri- 
bution to  the  total  attenuation  with  altitude. 

It  is  also  of  some  interest  to  have  information  concerning  both 
upgoing  and  downgoing  waves.  To  understand  how  it  is  achieved  a brief 
review  of  the  mechanics  of  the  full  wave  fields  program  will  first  be 
given. 

In  a horizontally  stratified  medium,  an  electromagnetic  wave  with 
time  and  x dependence  of  the  form  exp(i (ojt-kSx))  and  invariant  in  the  y 
direction,  has  an  altitude  (z)  dependence  determined  by  the  coupled  first 
order  differential  equations: 
de 

37  = (6) 

where  e is  a column  vector  of  field  components  given  by: 


it  • 


T is  a 4 X 4 matrix  which  depends  upon  the  incidence  angle,  o,  and  upon 
the  susceptibility  matrix  elements  (see  Clemmow  and  Heading^^  or  Budden^^ 

- Eq  18.7).  The  latter  of  course  depends  upon  z through  the  electron 
and  ion  densities  and  collision  frequencies.  They  also  depend  upon  both 
the  magnitude  and  orientation  of  the  geomagnetic  field. 

Various  computational  techniques  are  available  for  solving  Eqs.  (6). 
Smith' $10  program  is  based  upon  the  Inoue-HorowitzU  method  whereas 
Shellman's  variant  which  has  been  used  in  the  present  study  is  based  on  a 
a Runge-Kutta  integration,  but  otherwise  incorporates  the  same  scaling  and 


other  features  as  Smith's  program.  The  integration  is  begun  high  in  the 

ionosphere,  where  downcoming  waves  are  assumed  small.  The  two  upgoitig 

solutions  are  integrated  to  the  ground,  combined  to  achieve  the  correct 

polarization  at  the  ground  (except  for  unusual  circumstances  discussed 

in  section  4 and  Appendix  C the  polarization  at  the  ground  in  the  lower 

ELF  band  may  be  regarded  as  vertical)  and  re-combined  at  all  heights  to 

yield  correct  field  strengths.  Once  e (Eq.  6)  has  been  obtained  it  is 

possible  to  decompose  it  at  any  altitude  into  four  magneto-ionic  modes 

associated  with  a homogeneous  medium  having  the  properties  appropriate 

to  that  altitude.  Calling  the  magneto-ionic  solution  vectors,  e may 

be  written  as: 

4 

c = S (8) 

i=l 

where  6^  are  constants  detennined  by 

fjj  =S"*c  und  S['*Sj  = 6jj  where  6 = Kroncckcr  Delta  (9) 

Effecting  this  decomposition  it  is  then  possible  to  write  at  any  altitude: 

(10) 

where  u and  d designate  up  and  down  travelling  waves.  Substituting  the 
up  or  down  travelling  fields  into  (4)  we  define: 

ffU.d  =1  (1^) 

It  should  be  noted  that  in  general 

+ (12) 

Nevertheless,  at  high  enough  altitudes  plots  of  the  z caiiponent  of  (11) 
can  serve  as  a crude  measure  of  where  the  downcoming  wave  becomes  unim- 
portant. A somewhat  more  useful  index,  for  our  purposes,  of  where  the 
downcoming  wave  becomes  unimportant  is  pi'edicated  on  the  fact  that  at 
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high  altitudes  the  radio  wave  field,  except  in  a very  narrow  angular  band 
about  the  magnetic  equator  where  quasi -transverse  propagation  applies, 
is  almost  totally  in  the  outgoing  whistler  mode.  In  the  whistler  mode, 
the  wave  is  circularly  polarized  such  that  E2  = iE^  and  JCg  = i3f^. 

Thus,  measures  of  departure  of  the  total  full  wave  field  components 
from  the  outgoing  whistler  mode  are  given  by  the  following  quantities: 

(13) 

lE2-iEfl  (14) 

b.  Reflection. 

Westcott^^”^^  in  a series  of  papers  has  explored  ionospheric  reflec- 
tion processes  in  the  VLF  band  using  the  Darwin-Hartree^^’^^  microscopic 
theory  of  re-radiation.  In  this  theory,  the  re-radiation  of  the  electron 
due  to  its  acceleration  produced  by  the  total  field  is  of  primary  impor- 
tance and  leads  quite  naturally  to  a conception  of  how  each  level  of  a 
horizontally  stratified  medium  contributes  to  the  final  reflected  field. 
Using  the  full  wave  program  as  described  previously,  the  total  fields 
are  available  throughout  the  ionosphere  and  examination  of  reflection 
processes  using  the  microscopic  theory  is  easily  implemented.  The 
essential  formulas  are  summarized  below. 

The  most  important  result  of  the  mic»*oscopic  theory  for  our  purposes 
is  that  the  reflected  3f|  field  at  the  ground  produced  by  the  dipole 
layering  up  to  height  z is  given  by: 

ik  /*2  ' 

^2“  (Sp3+CP,)e'^‘‘^'^<ise  (15) 
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where  the  P-j's  are  polarization  vector  components  given  by  Eq.  (3).  By 
utilizing  the  relationship: 

(16) 


P 

Co 


V(y  • El-V*-  E -E 


the  integral  in  Eq.  (15)  may  be  carried  out  with  the  result: 
^2  = 7 [(-3f2  + E,/C)e“'''Czj^ 


(17) 


In  the  right  hand  side  of  Eq.  (17),  the3('s  and  E's  are  the  full  wave 
solutions  and  so  the  right  hand  side  is  readily  calculable  when  the 
latter  are  available.  Taking  the  derivative  of  Eq.  (15)  gives: 


dJf!, 

dz 


(SP3tCP,)c-iM>' 


(18) 


Equation  (18)  states  that  the  differential  field  re-radiated  by  a layer 
of  thickness  Az  at  height  z is  (d3(2/dz)A2  where  the  quantity  dH^/dz  is 
the  reflection  density  for  the  y component  of  the  magnetic  field.  Also, 
K2  given  by  Eq.  (17)  is  the  integrated  dipole  layer  contribution  to  the 
reflected  y component  of  the  magnetic  field. 
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A number  of  calculations  have  been  made  to  illustrate  the  type  of 
information  available  from  full  wave  outputs  as  regards  ionospheric 
absorption  and  reflection  features  of  ELF  waves.  Many  of  the  calculations 
have  been  made  using  the  ambient  day  and  night  ionospheric  models  schema- 
tized in  Fig.  1.  The  data  for  the  region  from  ground  to  100  km  are 
from  Knapp  and  those  above  100  km  are  from  Johnson.  Some  adjustments 
of  the  two  sets  of  data  .re  required  to  avoid  discontinuities  at  100  km. 

A mass  of  32  AMU  is  assigned  to  both  the  positive  and  negative  ions. 

The  collision  frequencies  as  a function  of  altitude  are  listed  in  Tables 
I and  II  for  daytime  and  nighttime,  respectively.  The  collision  freq- 
uencies vary  exponentially  with  altitude  between  the  tabulated  increments. 

No  claim  is  made  for  the  correctness  of  these  profiles.  Assignment 
of  a single  atomic  mass  unit  to  the  ion  population  is  clearly  a compro- 
mise and  the  electron  and  ion  densities  are  at  best,  only  reasonable 
representations  of  mean  conditions  taken  on  a global  basis.  It  is  known 
that  the  profiles  yield  attenuation  rates  and  field  predictions  in  reason- 
able agreement  with  the  Sanguine  and  Seafarer  test  measurements  but  at 
the  same  time  seemingly  suffer  the  shortcoming  of  overestimating  the 
attenuation  rate  in  the  neighboriiood  of  the  fundamental  Schumann  resonance. 

For  those  features  of  reflection  and  absorption  which  seem  to  be 
particularly  sensitive  to  n«)del  choice,  a number  of  other  models  will  be 
examined  as  well.  One  of  these  will  be  the  ambient  nighttime  model  of 
Fig.  1 with  the  sporadic  E layer  superimposed.  It  is  a layer  about  5 km 
thick  at  the  top  and  about  IS  km  thick  at  the  base  with  a peak  electron 
density  enhancement  factor  of  about  30  over  that  of  the  ancient.  The 


11 


Ftg.  1.  Day  and  Night  Ambient  Profiles  and  Nighttime  Sporadic 


TABLE  1. 


Daytime  Ionospheric  Electron  and  Ion  Collision 
Frequency  in  Seconds"^  Versus  Altitude 


Altitude 

(km) 

Electrons 

Positive 

Ions 

Negative 

Ions 

150. 

120. 

100. 

0.0 

i.eoxio^ 

i.ooxio; 

3.90x10^, 

4. 30x1011 

4.50x10} 

3.00x10^ 

8.00x103 

1.07xl0lO 

4.50xl0} 

3.00x10^ 

8.00x103 

1.07xl0l0 

TABLE  2. 


Nighttime  Ionospheric  Electron  and  Ion  Collision 
Frequency  in  Seconds”^  Versus  Altitude 


Altitude 

(km) 

Electrons 

Positive 

Ions 

Negative 

Ions 

250. 

1. 05x10^ 

4.50x10-1 

4.50x10-1 

9.00x10-1 

225. 

3. 50x101 
3.00x10; 

9.00x10-1 

220. 

1.00 

1.00 

210 

3.30xl0j 

4.50xl0l 

1.30 

1.30 

200. 

2.00  , 
4.50xl0j 

2.00  , 
4.50xl0l 

150. 

1.60x10? 

120. 

1.00x10} 

3.00x10; 

3.00x10; 

100. 

3.90x10}, 

8.00x10?- 

8.00x10?^ 

0.0 

4.30x1011 

1.07xl0l0 

1.07xl0l0 

layer  was  suggested  in  a private  communication  from  T.  R.  Larsen  of  the 
Norwegian  Defence  Research  Establishment.  During  the  course  of 
this  section,  description  of,  and  results  for  the  other  ionospheres 
will  be  given.  Although  finite  ground  conductivities  have  been  used  in 
the  calculations,  for  all  practical  purposes  the  results  presented  can 
be  thought  of  as  appropriate  for  an  infinitely  conducting  ground. 

Figures  2 through  4 show  full  wave  absorption  results  for  the  day- 
time ambient  profile  of  Fig.  1.  Figure  2 is  for  75  Hz,  Fig.  3 for  45  Hz 
and  Fig.  4 for  7.5  Hz  and  all  of  the  curves  are  for  an  azimuth  of 
propagation  of  270°  relative  to  magnetic  north,  a dip  of  60°  and  a 
geotiiagnetic  field  of  0.5  Gauss.  Height  in  kilometers  is  given  on  the 
vertical  axis,  relative  heating  loss  per  unit  volume  (Eq.  2)  on  the  lower 
horizontal  axis  and  integrated  absorption  measured  in  units  of  attenuation 
(i.e.  dB/1000  km  - Eq  5)  is  shown  on  the  upper  horizontal  axis.  Also 
shown  are  the  total  attenuations.  The  latter  are  certainly  compatible 
with  available  measurements.  Features  common  to  the  three  figures  are 
the  low  and  high  altitude  absorption  peaks.  The  low  altitude  absorption 
peaks  occur  at  “55  km  for  75  and  45  Hz  and  at  “51  km  for  7.5  Hz.  If 
the  profile  densities  and  collision  frequencies  vary  exponentially  in 
the  neighborhood  of  the  absorption  peaks  and  if  the  scale  heights  are 
invariant  over  the  line  profile  then  it  is  shown  in  Appendix  A that  tlHj 
absorption  line  shape  may  be  approximated  by  the  torentzian  line  shape 
given  by  Eq.  (A-8).  Also,  the  position  of  the  peak  is  given  approxiiMtely 
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HEIGHT  (KM) 
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Fig.  2.  Day  Absorption  and  Attenuation  at  75  hz.  Azimuth  = 270®,  Dip  = 
60®,  Geomagnetic  Field  = 0.5  Gauss, 
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Day  Absorption  and  Attenuation  at  45  Hz.  Azimuth  = 270 
Dip  = 60°,  Geomagnetic  Field  = 0.5  Gauss. 
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FiQo  4o  Day  Absorption  and  Attenuation  at  7.5  Hz.  Azimuth  = 270°, 
Dip  = 60°,  Geomagnetic  Field  = 0.5  Gauss. 
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by  Eq.  (A-6).  Comparisons  of  the  full  wave  low  altitude  absorption  lines 
(normalized  to  unity)  with  the  Lorentzian  lines  are  given  in  Figs.  5 
through  7 for  75,  45  and  7.5  Hz.  The  first  thing  to  observe  is  that  the 
altitude  where  the  low  altitude  absorption  peaks  is  quite  well  predicted 
by  Eq.  (A-6).  Apart  from  the  wings  of  the  lines,  the  agreement  between 
the  full  wave  output  and  the  Lorentzian  line  shape  given  by  Eq,  (A-8)  is 
excellent  for  the  75  Hz  and  45  Hz  cases  but  is  poor  for  the  7.5  Hz  case. 
This  poor  agreement  results  from  several  causes.  First,  the  half  width 
of  the  line  is  sufficiently  broad  that  the  assumption  of  invariant  scale 
heights  over  the  line  is  violated.  Second,  the  asymmetry  of  the  line 
could  only  be  accounted  for  by  retaining  higher  order  terms  in  the  ex- 
pansion of  the  denominators  of  Eqs,  (A-7)  and  (A-8). 

The  integrated  effect  of  the  low  altitude  absorption  line  (as  mea- 
sured in  units  of  attenuation)  has  also  been  approximated  in  Appendix  A 
by  Eq.  (A-13).  In  particular,  the  attenuation  will  be  seen  from  (A-13) 
to  be  particularly  sensitive  to  the  species  scale  heights  and  their 
collision  fj'equency  scale  heights  in  the  neighborhood  of  the  absorption 
peak.  From  Figures  Z through  4 It  will  be  seen  that  the  low  altitude 
absorption  contributes  about  0.4  dU/lOOO  km  to  the  attenuations  at  75  Hz, 
about  0.3  dB/1000  km  at  45  Hz  and  about  0.1  dB/1000  km  at  7,S  Hz. 

The  high  altitude  abso»i)tion  peaks  occur  at  about  114  kni  for  75  Hz, 
about  116  km  for  45  Hz  and  about  121  km  for  7.5  Hz.  If  the  profile  den- 
sities and  collision  frequencies  vary  exponentially  in  the  neighborhood 
of  the  absorption  peaks  and  if  the  scale  heights  ai*e  invariant  over  the 
line  profile  then  it  is  shown  in  Appendix  0 that  the  absoriition  line 
shape  may  be  approximated  by  the  Gaussian  line  shape  given  by  Eq.  (6-13). 
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Fig,  6.  Low  Altitude  Daytime  Absorption  Comparisons  at  45  Hz, 
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Low  Altitude  Oaytinie  Absorption  Comparisons  at  7o5  Hz. 


Also,  the  position  of  the  peak  is  given  approximately  by  Eq.  (B-12). 
Comparisons  of  the  full  wave  absorption  lines  (normalized  to  unity)  with 
the  Gaussian  lines  are  given  in  Figure  8 through  10  for  75,  45  and 
7.5  Hz.  Observe  that  the  altitude  where  the  high  altitude  absorption 
peaks  is  quite  well  predicted  by  Eq.  (B-12).  Apart  from  the  wings  of 
the  lines,  the  agreement  between  the  full  wave  output  and  the  Gaussian 
line  shape  given  by  Eq.  (B-13)  is  good  for  the  75  and  45  Hz  cases  but 
is  poor  for  the  7.5  :iz  case.  As  was  the  case  with  the  low  altitude 
absorption  line  comparisons,  this  poor  agreement  results  from  the  fact 
that  the  line  is  sufficiently  broad  to  invalidate  the  assumption  of 
invariant  scale  heights.  Furthermore,  the  marked  asymmetry  of  the  7.5 
Hz  high  altitude  line  could  only  be  accounted  for  by  retaining  higher 
order  teniis  in  the  expansion  of  the  exponent  of  Eq.  (B-13). 

From  Figures  2 through  4 it  will  be  seen  that  the  high  altitude 
absorption  contributes  about  0.2  dB/1000  km  to  the  attenuation  at  75  Hz, 
about  0.4  dB/1000  km  at  45  Hz  and  about  0,06  d8/1000  km  at  7.5  Hz  (of 
course  it  may  be  one's  preference  to  include  these  losses  as  part  of  a 
leakage  loss  out  of  the  top  of  the  guide).  Unlike  the  low  altitude 
absorption,  it  is  not  an  easy  task  to  give  an  approximate  estimate  of 
the  integrated  absorption  associated  with  the  high  altitude  line  because 
it  depends  upon  the  transmission  through  the  D and  E layers.  Similarly, 
the  intermediate  altitude  absorption  lines  cojianon  to  Figures  2 through  4 
are  not  readily  amenable  to  analysis  since  they  occur  in  a r’egion  over 
which  upgoiruj  and  d(Mngoing  magneto- ionic  components  coexist  in  variable 
proportions. 

Figures  11  through  13  show  full  wave  absorption  results  for  the 
nighttime  ambient  profile  of  Figure  1.  Figure  U is  for  75  Hz,  Figure  12 
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Fig.  9.  High  Altitude  Daytime  Absorption  Comparisons  at  45  Hz 
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Fig.  11.  Night  Absorption  and  Attenuation  - 75  Hz.  Azimuth  - 270®, 
Dip  » 60°.  Geomagnetic  Field  * 0.5  Gauss. 


26 


HEIGHT 


Fig.  12.  Night  Absorption  and  Attenuation  - 45  Hz,  Azimuth  = 270' 
Dip  ® 60® » Geomagnetic  Field  *»  0«5  Gauss. 
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FI9.  13.  Night  Absorption  and  Attenuation  - 7.5  Hz.  Aziaiuth  2 0 . 
Dip  =>  60“,  Geomagnetic  Field  “ 0,5  Gauss, 
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for  45  Hz  and  Figure  13  for  7.5  Hz  and  all  of  the  curves  are  for  an 
azimuth  of  propagation  of  270°  relative  to  magnetic  north  a dip 
of  60°  and  a geomagnetic  field  strength  of  0.5  Gauss.  The  curves  show 
the  same  quantities  as  those  of  Figures  2 through  4 for  the  daytime 
profiles.  The  total  attenuations  shown  on  the  figures  are  in  the  45  and 
75  Hz  cases,  compatible  with  available  measurements  whereas  the  attenua- 
tion rate  of  1.3  dB/1000  km  for  7.5  Hz  seems  quite  inconsistent  with  the 
observation  of  Schumann  resonances,  although  the  latter,  of  course,  de- 
pend upon  a global  fverage  for  which  the  mid  latitude  nighttime  proto- 
type of  Figure  1 may  play  a relatively  minor  role.  At  any  rate  the 
nighttime  profile  of  Fiyiue  1 can  at  least  be  thought  of  as  empirical 
for  Seafarer  band  pruuictions. 

As  was  the  case  for  the  daytime  profile.  Figures  11  through  13 
show  the  common  features  of  low  and  high  altitude  absorptions.  The  low 
altitude  absorption  peaks  occur  at  “72  km  for  75  and  45  Hz  and  at  “62  km 
for  7.5  Hz.  As  discussed  in  Appendix  A,  the  analysis  given  there,  with 
very  minor  alteration,  can  also  be  applied  at  least  approximately  to  the 
nighttime  low  altitude  absorption.  Comparisons  of  the  full  wave  low 
altitude  absorption  line  locations  and  shapes  (nonnalized  to  unity)  with 
the  Lorentzian  line  shapes  of  Appendix  A are  given  in  Figures  14  through 
16.  As  was  the  case  with  the  daytime  results,  the  locations  of  the  peaks 
are  quite  well  predicted  by  the  simple  theory  of  Appendix  A.  However, 
the  full  wave  low  altitude  line  shapes  for  the  nighttime  case  study  show 
more  asymmetry,  particularly  at  76  and  45  Hz,  than  did  the  daytinse  results. 
Still,  If  wing  effects  are  discounted  the  agreement  at  all  frequencies 
is  fairly  good,  from  Figures  11  through  13  it  will  be  seen  that  the 
contribution  of  the  low  altitude  absorption  line  is  “0.5  dl3/1000  km  at 
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Fig.  15.  Lom  Altitude  Nighttime  Absorption  Comparisons  at  45  Hz. 


75  Hz,  =0.4  dB/1000  km  at  45  Hz  and  =0.4  dB/1000  km  at  7.5  Hz. 

The  high  altitude  absorption  peaks  occur  at  =113  km  for  75  Hz,  at 
=115  km  for  45  Hz  and  at  =119  km  for  7.5  Hz,  Comparisons  of  the  full 
wave  line  shapes  (normalized  to  unity)  with  the  Gaussian  line  shapes  pre- 
dicted by  the  simple  theory  of  Appendix  B are  shown  in  Figures  17,  18 
and  19.  Again  if  wing  effects  are  discounted,  the  agreement  is  quite 
good  with  the  7.5  Hz  being  again  the  worst  of  the  three  frequencies 
considered. 

From  Figures  11  through  13  it  will  be  seen  that  the  contribution 
of  the  high  altitude  absorption  to  the  attenuation  is  =0.1  dB/1000  km 
at  75  and  45  Hz  and  =0.4  dB/1000  km  at  7.5  Hz.  The  latter  figure  is  an 
indication  of  the  quite  appreciable  energy  penetrating  the  nighttime 
ambient  ionosphere  of  Figure  1 to  altitudes  of  120  km  and  above.  Further 
evidence  of  this  will  be  forthcoming. 

Figures  20  and  21  show  results  relevant  to  reflection  properties  of 
the  nighttime  ambient  profiles  of  Figure  1 as  looked  at  from  the  point 
of  view  of  the  Darwin-Hartree^^’^^  re-radiation  theory  as  discussed  in 
section  3.  The  figures  apply  to  the  same  geomagnetic  conditions  as 
Figures  11  through  13  and  more  particularly  Figure  20  is  for  75  Hz  and 
Figure  21  is  for  45  Hz.  The  vertical  scale  in  each  instance  is  the  height 
in  km.  The  upper  horizontal  scale  applies  to  the  solid  curve  (Eq.  (17) 
of  section  3)  and  to  the  dotted  curve  (£q.  (13)  of  section  3).  The  lower 
horizontal  scale  applies  to  the  dashed  curve  (Eq.  (4)  of  section  3). 

The  solid  curve  represents  the  buildup  of  the  magnitude  of  the  y compon- 
ent of  the  reflected  3C  field  at  the  ground  due  to  the  integrated  effect 
of  the  dipole  layer  structure  below  the  height  indicated  on  the  vertical 
axis.  Because  of  the  noniialization  used,  the  value  of  the  magnitude  of 
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Fig,  17.  High  Altitude  Nighttime  Absorption  Comparisons  at  75  Hz 
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Fig,  19,  High  Altitude  Nighttime  Absorption  Comparisons  at  7,5  Hz 


the  y component  of  the  reflected  Jf  field  at  the  ground  associated  with 
the  re-radiation  from  the  total  ionosphere  is  0.5.  The  dotted  curve  is 
the  magnitude  of  the  dashed  curve  is  the  magnitude  of  £2“^^! 

where  the  field  components  are  either  full  wave  components  or  decompos- 
itions thereof  as  discussed  in  section  3.  Thus  the  dotted  and  dashed 
curves  may  be  thought  of  as  measures  of  the  departure  of  the  full  wave 
solutions  from  the  outgoing  whistler  mode,  measures  which  are  meaningful 
principally  at  high  altitudes,  where  the  total  field  is  expected  to  be 
in  the  outgoing  whistler  mode. 

The  solid  curve  of  Figure  (20)  illustrates,  then,  that  the  buildup 
of  the  magnitude  of  the  y component  of  the  reflected  field  at  the  ground 
begins  slowly  around  70  km,  is  quite  rapid  in  the  95  to  120  km  range  and 
attains  in  magnitude  its  fully  developed  value  of  0.5  at  about  133  km. 

A measure  of  the  integrated  contribution  of  the  ionosphere  above  this 
height  to  the  y component  of  the  reflected  iJf  field  at  the  ground  is  the 
departure  of  its  phase  from  its  fully  developed  value  of  0.0°.  This 
departure  in  the  case  of  Figure  (20)is  known  from  the  full  wave  output 
to  be  only  about  11°  and  so  we  interpret  this  as  indicating  that  the 
reflection  as  looked  at  from  the  re-radiation  point  of  view  is  in  large 
measure  developed  by  the  ionosphere  below  the  133  km  height,  despite  the 
fact  that  the  solid  curve  does  not  asymptotically  level  off  beyond  there 
to  its  final  value  of  0.5.  A necessary  (but  not  sufficient)  condition 
for  this  to  be  a valid  interpretation  is  that  there  be  no  appreciable 
downcoming  wave  above  133  km.  That  this  is  in  fact  the  case  is  indicated 
by  the  dashed  and  dotted  curves.  These  curves  show  that  the  fields 
above  about  130  km  are  very  nearly  described  by  the  outgoing  whistler 
mode.  The  large  amplitude  oscillatory  phase  of  the  solid  curve  above 
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200  km  is  due  to  this  whistler  mode  penetration  coupled  with  the  large 
increase  in  the  electron  density  associated  with  the  F layer  (see  Figure  1). 

The  interpretation  of  Figure  (21)  which  applies  to  45  Hz,  not  surpris- 
ingly, is  very  much  the  same  as  that  of  Figure  (20).  In  this  case  the 
magnitude  of  the  y component  of  the  reflected  3f  field  at  the  ground 
attains  the  value  of  0.5  at  about  146  km.  At  that  altitude  the  phase 
departure  of  the  y component  of  the  reflected  3f  field  at  the  ground  is 
about  10*^  from  its  final  value  of  0.0°.  Thus,  as  before,  it  is  concluded 
that  the  reflection  from  the  re-radiation  point  of  view  is  in  large 
measure  developed  by  the  ionosphere  below  about  146  km. 

An  alternative  display  of  the  levels  contributing  to  the  reflection 
is  given  in  Figures  (22)  and  (23)  where  the  magnitude  of  the 363  t^eflected 
field  density  dJfpdz  as  defined  by  Eq.  (18)  has  been  plotted.  Also  shown 
on  the  plots  are  the  relative  heating  losses  previously  shown  in  Figures 
11  and  12.  The  large  amplitude  increase  inldUf^/dzl  above  200  km  is  dis- 
counted for  the  reasons  cited  above.  It  is  due  to  polarizations  assoc- 
iated with  outgoing  (WKB  like)  wave  penetrating  the  F region  and  on  the 
basis  of  the  preceding  discussion  is  known  to  contribute  little  to  the 
integrated  reflected  3(3  ^ield  at  the  ground.  Thus  the  major  reflection 
comes  in  the  case  of  the  75  Hz  wave  from  the  region  about  the  peak  at 
about  103  km.  In  the  case  of  the  45  Hz  wave  the  major  reflection  comes  from 
the  region  about  the  peak  at  about  102  km.  The  half  widths  in  both  cases 
are  about  20  km.  Figures  (22)  and  (23)  illustrate  quite  graphically  the 
significant  regions  as  regards  both  reflection  and  absorption. 

Figure  (24)  shows  reflection  infonnation  for  the  nighttime  profile 
at  7.5  Hz.  The  solid  curve  illustrates  the  buildup  of  the  magnitude  of 
the  y canponent  of  the  ♦’eflected  3(  field  at  the  ground  due  to  the 
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integrated  effect  of  the  dipole  layer  structure  below  the  height  indicated 
on  the  vertical  axis.  Because  the  modal  height  gain  at  the  starting 
height  of  250  km  is  in  this  instance  composed  of  two  outgoing  magneto- 
ionic components,  it  is  not  convenient  to  present  field  results  as  for 
the  75  and  45  Hz  cases  in  terms  of  departures  from  the  outgoing  whistler 
mode.  Rather,  we  present  in  Figure  (24)  for  altitudes  above  100  km  the 
upgoing  and  downgoing  Poynting  flux  as  defined  by  Eqs.  (11).  These  curves 
show  that  the  upgoing  and  downgoing  fluxes  are  of  the  same  order  of 
magnitude  throughout  much  of  this  height  range,  thereby  indicating  that 
the  final  values  depend  in  some  measure  on  reflections  occurring  at  the 
uppermost  altitudes.  The  solid  curve  solidifies  this  conclusion  in  the 
following  sense.  The  solid  curve  achieves  the  value  of  0.5  at  an  alti- 
tude of  about  115  km,  however  the  phase  of  the  y component  of  the  re- 
flected field  at  the  ground  departs  from  its  final  value  of  0.0°  by 
about  113°.  We  conclude  then  that  the  integrated  re-radiation  above 
this  level  is  significant.  Apparently  because  of  interference  between 
reflections  from  disparate  ionospheric  regions  of  the  ionosphere  (see 
e.g.,  the  steep  gradients  in  the  solid  curve  of  Figure  (24)  around  110 
and  240  km)  the  normal  mode  polarization  is  not  TM  but  has  in  fact  a 
substantial  TE  admixture.  A measure  of  this  mixing  is  the  magnitude  of 
the  ratio  of  3Cj/3(2  at  the  ground.  This  value  in  the  present  instance 
is  about  1.3. 

To  understand  how  interference  frojit  ionospheric  reflecting  layers 
can  lead  to  polarization  mixing  we  will  turn  now  to  a discussion  of  the 
ambient  nighttime  profile  of  Figure  (1)  with  the  superimposed  sporadic  E 
layer.  Recall  that  this  Is  a layer  around  120  km,  about  5 km  thick  at 


the  top  and  about  15  km  thick  at  the  base  with  a peak  electron  density 
enhancement  factor  of  about  30  over  that  of  the  nighttime  ambient. 

Figure  (25)  shows  the  attenuation  versus  frequency  for  both  the  night- 
time ambient  and  for  the  nighttime  ambient  with  the  superimposed  sporadic 
E.  As  will  be  seen,  there  exists  a strong  resonance  which  at  its  peak 
represents  about  an  order  of  magnitude  enchancement  over  the  ambient 
attenuation.  Also  it  will  be  seen  that  the  attenuation  peaks  at  about 
78  Hz.  Figure  (26)  shows  the  relative  heating  loss  per  unit  volume  and 
its  integrated  effect  in  attenuation  units  for  this  case  at  75  Hz.  The 
curve  shows  the  low  altitude  absorption  at  about  72  km,  a much  amplified 
(over  the  nighttime  ambient  profile  alone  - see  Figure  (11))  high  alti- 
tude absorption  at  about  105  km  and  an  absorption  associated  with  the 
sporadic  E layer.  The  low  altitude  absorption  contributes  less  than  1 
dB/1000  km  to  the  attenuation.  The  high  altitude  absorption  contributes 
about  4 dB/1000  km  and  the  sporadic  E layer  itself  contributes  about  2.5 
dB/1000  km.  The  remainder  of  about  2.5  dB/1000  km  can  be  attributed  to 
leakage  beyond  the  sporadic  E layer. 

Figure  (27)  shows  again  the  relative  heating  loss  but  along  with 
the  magnitude  of  the  reflectedJfg  field  density  defined  by  Eq.  (18)  of 
section  3.  The  increase  in  the  latter  curve  in  the  neighborhood  of  240 
km  is  discounted  because  the  rf  field  above  the  sporadic  E layer  is  prin- 
cipally in  the  outgoing  whistler  mode  (See  Figure  (28)  where  the  quanti- 
ties IWg-lwVl  and  lEg-Itjl  have  been  plotted).  Thus  Figure  (27)  shows  a 
reflection  region  centered  about  100  km  and  another  centered  at  the 
sporadic  E layer.  To  better  understand  the  origin  of  the  resonance  be- 
havior of  Figure  (26),  the  preceding  tMo  layer  reflection  behavior 
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suggests  the  simplified  model  shown  in  Figure  (29).  It  consists  of 
vacuum  below  level  2,  equal  and  homogeneous  electron  and  ion  densities 
N2»  N3  and  in  regions  2,  3,  and  4 respectively.  The  electron  and  ion 
collision  frequencies  are  also  assumed  constant  in  the  different  regions. 
The  model  also  assumes  that  quasi  longitudinal  propagation  applies  in  the 
altitude  region  above  z.  The  analysis  for  this  model  is  given  in  Appen- 
dix C.  For  parameters  like  those  appropriate  to  Figure  (27),  the  results 
of  Appendix  C suggest  that  a primary  condition  for  an  accentuated  reson- 
ance is  that  the  region  between  the  two  reflecting  levels  (i.e.,  between 
and  Z2  in  Figure  (29))  behave  like  a quarter  wave  plate  and  that  simul- 
taneously the  sporadic  E layer  also  have  a thickness  equal  to  one  quarter 
integer  multiples  of  the  local  wavelength,  thereby  reducing  leakage 
through  the  layer.  The  resonance  is  shifted  in  frequency  and  consider- 
ably diminished  if  the  thickness  of  the  sporadic  E layer  is  half  integer 
multiples  of  the  local  wavelength,  a condition  conducive  to  leakage  out 
of  the  top  of  the  layer.  Location  and  strength  of  the  resonance  varies 
between  the  above  limits  as  the  thickness  of  the  sporadic  E layer  changes 
* from  quarter  wave  length  behavior  to  half  wave  length  behavior.  These 

1 features  are  illustrated  in  Figure  (30)  which  shows  attenuation  versus 

I frequency  for  the  following  parameters  of  Figure  (29)  (subscripts  e and 

I i stand  for  electrons  and  ions  respectively): 

I Ng  = N4  = 2 X 10^/cn»^  . N3  = 6 X loVcm^ 

f'  d 1 

I ^e2  = ^*e3  ^e4 
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Z2~^l  ~ 20  km  Quarter  wavelength  condition  at  123  Hz 
23-Z2  = 3.43  km  Quarter  wavelength  condition  at  123  Hz 
Z3-Z2  = 5.71  km 

Z3-Z2  = 6.86  km  Half  wavelength  condition  at  123  Hz 
The  heavy  curves  in  Figure  (30)  are  full  wave  output  and  the  light  curves 
are  model  results  calculated  using  Eqs.  C.IO,  C.23  and  C.25  of  Appendix 
C.  It  will  be  seen  that  the  full  wave  results  and  the  simplified  model 
results  are  in  excellent  agreement.  It  will  further  be  noted  that  the 
resonance  is  shifted  and  diminished,  as  stated  above,  when  region  3 has 
a thickness  equal  to  one  half  of  the  local  wavelength.  Figure  (31) 
shows  the  relative  heating  loss  per  unit  volume  and  its  integrated 
effect  in  attenuation  units  at  123  Hz  for  the  case  when  region  3 behaves 
like  a quarter  wave  plate.  The  absorption  below  120  km  contributes  close 
to  5 dB/lOOO  km  to  the  attenuation  and  region  3 contributes  by  virtue  of 
absorption  about  1 dB/100  km.  The  remainder  of  about  3 dB/1000  km  can  be 
attributed  to  leakage  into  region  4.  These  numbers  are  similar  to  those 
quoted  in  connection  with  Figure  (26)  and  in  fact,  not  surprisingly,  a 
striking  similarity  will  be  noted  between  Figures  (26)  and  (31).  Figure 
(32)  shows  the  relative  heating  loss  per  unit  volume  and  its  integ'*ated 
effect  in  attenuation  units  at  SO  Kz  for  the  case  when  the  thickness  of 
region  3 is  one  half  of  the  local  wavelength.  In  this  case  the  reflection 
by  the  sporadic  E layer  into  the  earth- ionosphere  waveguide  is  much  reduced 
from  the  previous  case  so  that  the  field  buildup  and  the  attenuation  in  the 
100  to  120  km  range  is  very  much  reduced  from  the  previous  case. 

Under  the  condition  of  strong  attenuation  resonance,  there  is  appreciable 
TE  admixture  In  the  normal  mode.  A measure  of  this  is  given  by  the  ratio  of 
at  the  ground.  This  ratio  is  for  the  conditions  appropriate  to  Figure 

(31). 
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= 0.747/-15.7°  from  full  wave  calculation 

7fj/3(’2  = 0.723/-12.5°  from  C.23  and  C.38 

Although  the  profiles  of  Figure  1 yield  modal  parameters  which  agree 
quite  well  with  available  data  in  the  Seafarer  band,  they  dc  seem  to  yield 
an  unreasonably  high  attenuation  rate  in  the  range  of  the  fundamental 
Schumann  frequency.  This  is  because  of  excessive  leakage  above  the  E 
layer.  Booker  and  Lefeuvre,  for  this  reason,  have  suggested  that  pro- 
files with  stronger  F.  region  ledges  may  be  more  reasonable  and  one  might 
legitimately  ask  if  the  resonance  behavior  shown  in  Figure  25  would  then 
be  wiped  out.  Comparison  of  such  a profile  (referred  to  as  the  B-L  pro- 
file) with  the  nighttime  profile  of  Figure  1 is  shown  in  Figure  (33), 

The  much  stronger  E region  ledge  of  the  B-L  profile  is  evident.  Figure 
34  shows  attenuation  rate  versus  frequency  for  the  nighttime  ambient 
profile  with  the  sporadic  E layer  superimposed  along  with  the  results  for 
the  B-L  profile  with  the  sporadic  E layer  and  still  another  curve  showing 
results  for  the  B-L  profile  with  the  sporadic  E layer  lowered  by  10  km 
(i.e.,  the  layer  is  at  about  110  km  rather  than  120  km).  It  will  be 
seen  that  the  resonance  is  most  pronounced  for  the  night  ambient  profile 
with  the  superimposed  sporadic  E layer.  The  reason  for  this  is,  as  just 
discussed,  that  the  quarter  wavelength  conditions  are  satisfied.  The 
resonance  is  diminished  and  shifted  out  of  the  Seafarer  band  for  the  case 
of  the  sporadic  E layer  at  120  km  superimposed  on  the  B-L  profile  but  is 
again  shifted  into  the  Seafarer  band  if  the  sporadic  E layer  is  lowered 
10  km  and  superimposed  on  the  B-L  profile.  The  significance  of  these 
results  Is  of  course  the  degree  of  sensitivity  to  the  ledge  ionisations. 

This  is  further  evidenced  by  Figure  35  which  again  shows  attenuation  rate 
versus  frequency.  In  this  case  the  sporadic  E layer  has  been  modeled  as 
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and  B-L  + Sporadic  E.  Azimuth  = 105°.  Dip  = 77°,  Geomagnetic 
Field  = 0.41  Gauss. 


a sharply  bounded  homogeneous  slab  starting  at  110  km  and  having  thick- 
ness of  2o5,  5,0,  and  10,0  km.  Again  the  effect  is  significant  in  the 
Seafarer  band  and  as  expected  from  the  previous  discussion  of  the  simpli- 
fied model  of  Figure  29  shows  marked  sensitivity  to  the  modeled  sporadic 
E layer  thickness. 

We  now  turn  to  a discussion  of  reflection  properties  associated 
with  the  daytime  profile  of  Figure  1.  In  particular  Figures  36  through 
38  show  results  relevant  to  the  re-radiation  theory  as  discussed  in 
section  3.  The  figures  apply  to  the  same  geomagnetic  conditions  as 
Figures  2 through  4 and  more  particularly  Figure  36  is  for  75  Hz,  Figure 
37  for  45  Hz  and  Figure  38  for  7.5  Hz.  The  upper  horizontal  scale  applies 
to  the  solid  curve  (Eq.  (17)  of  section  3)  and  to  the  dotted  curve  (Eq.. 
(13)  of  section  3).  The  lower  horizontal  scale  applies  to  the  dashed 
curve  (Eq.  (14)  of  section  3).  The  solid  curve  represents  the  buildup 
of  the  magnitude  of  the  reflected  y cojiiponent  of  the  iJf  field  at  the 
ground  due  to  the  integrated  effect  of  the  dipole  layer  structu>'e  below 
the  height  indicated  on  the  vertical  axis.  Because  of  the  nonnalization 
used  the  value  of  the  magnitude  of  Mp  associated  with  the  re-radiation 
from  the  total  ionosphere  is  0.5.  The  dotted  curve  is  the  magnitude  of 
3fp-lK^and  the  dashed  curve  is  the  magnitude  of  E2-iEj  where  the  field 
components  ar“e  either  full  wave  components  or  decompositions  the>'eof 
as  discussed  in  section  3. 

The  solid  curve  of  Fig.  36  illustrates,  then,  that  the  buildup  of 
tlH?  magnitude  of  the  y cwnponent  of  the  reflected  field  at  the  ground 
begins  slowly  around  55  km  and  achieves  its  final  value  in  magnitude  of 
0.5  at  about  90  km.  Since  the  total  field  above  about  95  km  is  almost 
totally  in  the  outgoing  whistler  mode,  the  large  amplitude  oscillations 
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Fig.  36.  Reflection  and  Field  Penetration  Properties  for  Day  at  75  Hz 
Azimuth  = 270“,  Dip  = 60“ , Geomagnetic  Field  = 0.5  Gauss, 
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in  the  solid  curve  above  this  height  are  due  almost  exclusively  to  the 
outgoing  whistler  mode  leakage  and  the  resulting  large  polarization 
associated  with  the  daytime  ionosphere.  On  the  basis  of  the  combined 
information  conveyed  by  the  solid,  dotted  and  dashed  curves  of  Fig.  36, 
the  conclusion  is  that  the  reflection  is  in  large  measure  accomplished 
by  the  dipole  structure  below  the  95  km  level.  Figure  37  illustrates 
the  same  features  as  those  just  discussed  ’n  connection  with  Fig.  36. 

Figure  38  shows  the  considerably  deeper  field  penetration  which  occurs 
for  7.5  Hz.  In  this  case  the  reflection  appears  to  be  determined  prin- 
cipally by  the  dipole  structure  below  about  120  km.  Figures  36  through 
38  all  indicate  a reflection  behavior  similar  to  that  associated  with 
incidence  on  a pure  dielectric  and  the  fact  that  FigSo  36  through  38  do 
not  indicate  clearcut  reflecting  levels  result  from  this  behavior.  The 
penetration  depths  of  95  km  and  120  km  quoted  above  can  best  be  likened 
to  the  notion  of  a skin  depth  although  the  ionosphere  is  behaving  in  some 
respects  more  like  a dielectric  reflector  rather  than  a metallic  reflector. 
There  are  circumstances  when  the  ionosphere  behaves  l^ke  a metallic 
reflector.  Absorption  properties  for  these  special  cases,  especially 
as  regards  the  influence  of  low  level  ions,  has  been  discussed  rather 
extensively  by  Field. One  of  the  special  cases  Is  when  propagation  in 
the  ionosphere  can  be  approximated  by  quasi  transverse  propagation.  This 
requirest  as  discussed  by  Booker  and  Lefuvi'e,  propagatio*i  in  a nar- 
row angular  sector  about  the  magnetic  equator  and  Figs.  39  through  41  show 
what  can  be  expected  under  those  circumstances.  The  figures  apply  to  the 
daytime  profile  under  geomagnetic  conditions  of  east  to  west  propagatioit 
at  the  magnetic  equator.  The  upper  scale  applies  to  the  solid  curve  and 
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Fig,  41 o Reflection  Density  and  Absorption  for  Day  at  the  Magnetic 
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Field  = 0.5  Gauss. 
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the  lower  scale  applies  to  the  dotted  curve.  The  solid  curve  is  the  rel- 
ative heating  loss  per  unit  volume  and  the  dotted  curve  is  the  magnitude 
of  the  ("eflected  field  density  (dJf^/dz)  and  shows  in  these  instances 
well  defined  reflection  levels  in  the  neighborhood  of  84  km  for  the  75 
and  45  Hz  cases  and  in  the  neighborhood  of  90  km  for  the  7.5  Hz  case. 

The  integrated  buildup  of  the  magnitude  of  the  y component  of  the  re- 
flected 3{'2  i'ield  at  the  ground  due  to  the  dipole  structure  beneath  the 
ordinate  value  (Eq.  (17)  of  section  3)  for  these  cases  is  shown  in  Figs. 

42  through  44.  The  leveling  off  of  the  curve  at  a value  of  0.5  should 
be  compared  with  the  oscillating  behavior  shown  in  Figs.  36  through  38. 

When  the  ionosphere  is  highly  depressed  it  again  acts  like  a met- 
tallic  reflector.  This  is  illustrated  in  Figs.  46  through  48.  These 
figures  are  for  the  hypothetically  depressed  ionosphere  shown  in  Fig,  45, 
The  curves  show  the  same  information  and  have  the  same  scale  legends  as 
Figs,  39  through  41,  They  show  clearcut  reflection  levels  at  about  65  Ion 
for  the  75  and  45  Hz  cases  and  a considerably  broader  reflection  stratum 
centered  about  70  km  for  the  7.5  Hz  case.  Note  also  that  in  these  cases 
there  are  two  relatively  low  altitude  absorption  lines,  the  lowest  being 
centered  about  40  km  (due  to  ions)  and  the  other  about  60  km.  Figs.  49 
through  51  show  the  integrated  buildup  of  the  magnitude  of  the  y conponent 
of  the  reflected  }(  field  at  the  ground  due  to  the  integrated  dipole  struc- 
ture of  the  ionosphere  below  the  ordinate  value.  Note  again  the  rather 
rapid  (cwnpared  to  Figs,  36  through  38)  asymptotic  approach  of  the  curves 
to  the  0,5  value. 
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Fig.  42o  Reflection  Buildup  for  Day  at  the  Magnetic  Equator  at  75  Hz 
Azimuth  = 270°,  Dip  = 0°,  Geomagnetic  Field  =0.5  Gauss, 
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rig,  47c  Reflection  Density  for  Depressed  Ionosphere  at  45  Hz, 
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Fig.  48o  Reflection  Density  for  Depressed  Ionosphere  at  7o5  HZo 
Azimuth  = 45°,  Dip  = 75°,  Geomagnetic  Field  = 0o5  Gauss 
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Fig.  49o  Reflection  Buildup  for  Depressed  Ionosphere  at  75  Hz, 

Azimuth  = 45®,  Dip  = 75°,  Geomagnetic  Field  = 0,5  Gauss, 
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5.  CONCLUSIONS 


Under  conditions  of  mid  latitude  propagation  the  full  wave  solutions 
for  the  day/night  ambient  profiles  of  Fig.  1 exhibit  a low  altitude  ab- 
sorption line,  due  principally  to  electrons,  which  peaks  at  about  55  km 
during  the  day  and  around  70  km  during  the  night  as  well  as  a high  alti- 
tude absorption  line,  due  to  ions,  which  peaks  around  115  km  for  both 
day  and  night.  Theoretical  analysis,  based  on  the  assumption  of  scale 
height  invariance  of  the  profile  species  and  their  collision  frequencies 
over  the  absorption  line,  suggests  that  the  low  altitude  line  shape  should 
be  essentially  Lorentzian,  whereas  the  high  altitude  absorption  line  shape 
should  be  approximately  Gaussian.  For  the  case  studies  examined,  agree- 
ment between  these  theoretical  line  shapes  and  the  fullwave  line  shapes 
is  reasonably  good  in  the  Seafarer  band  (see  Figs.  5,  6,  8,  9,  14,  15,  17, 
18)  but  generally  poor  at  the  fundamental  Schumann  frequency  because  of 
the  large  halfwidths  (see  Figs.  7,  lO,  16,  19). 

Considerable  attention  has  been  given  to  the  question  of  resonances 
in  the  attenuation  associated  with  nighttime  sporadic  E layering.  Fig. 

25  shows  that  the  sporadic  E layer  of  Fig.  (1)  causes  an  order  of  magni- 
tude enhancement  in  the  attenuation  at  78  Hz  over  the  nighttime  ambient 
f'esult.  About  40^  of  the  enhanced  value  is  due  to  the  high  altitude 
absorption,  about  25%  to  absorption  which  occurs  in  the  layer  itself, 
about  25%  to  leakage  above  the  sporadic  E layer  and  about  10%  to  the  low 
altitude  absorption  (see  Fig.  26).  It  is  concluded  from  the  dipole  re- 
radiation theory  that  most  of  the  reflection  from  the  nighttime  ambient 
profile  close  to  the  resonance  frequency  comes  from  a region  about  20  km 
thick  centered  at  100  km.  The  sporadic  E layer  provides  a strong 
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reflection  level  at  120  km.  This  two  level  reflection  characteristic 
suggested  the  simplified  model  of  Fig.  (29)  and  on  the  basis  of  this 
simplified  model  it  has  been  concluded  that  the  fundamental  attenuation 
resonance  is  most  accentuated  when  the  separation  between  the  two  re- 
flecting layers  is  one  quarter  of  the  local  wavelength  and  when,  simul- 
taneously, the  sporadic  E layer  thickness  is  also  a quarter  of  the  local 
wavelength.  Both  of  these  conditions  are  conducive  to  waveguiding  be- 
tween the  two  reflection  levels.  Although  the  effect  is  diminished  for 

an  ambient  nighttime  profile  with  a stronger  E region  ledge  (as  suggested 

23 

by  Booker  and  Lefeuvre  ) it  appears  that  the  resonance  could  still  be 
of  practical  significance  (see  Figs,  34,  35), 

Subject  to  conditions  of  strong  resonance,  unlike  the  normal  situ- 
ation, the  normal  mode  polarization  contains  significant  amount  of  TE 
admixture.  Thus,  in  future  measurements  of  anomalous  propagation  of  the 
type  discussed  in  references  8 and  9 it  is  strongly  recommended  that  a 
provision  be  made  for  measuring  the  extent  of  TE  contaminant. 

In  the  Seafarer  band  field  penetration  and  re-radiation  for  the 
day  ambient  profile  of  Fig.  1 at  mid  latitudes,  occurs  up  to  at  least 
altitudes  of  90  kin  (see  Figs.  36,  37).  For  the  nighttime  ambient  profile 
of  Fig.  1 at  mid  latitudes,  significant  field  penetration  and  re-radiation 
occurs  to  altitudes  of  about  130  km  (Fig.  20)  for  75  Hz  and  to  altitudes 
of  about  145  km  (Fig.  21)  at  45  Hz. 

For  the  daytime  profile  the  re-radiation  theory  does  not  give  a 
clearcut  picture  of  the  level  where  maximum  re-radiation  occurs,  but 
rather  yields  an  indecisive  picture  somewhat  reminiscent  of  incidenpe 
on  a pure  dielectric.  Unlike  the  daytime  case,  the  level  of  maximum 
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reflection  based  on  the  re-radiation  theory  for  the  satellite  nighttime 
profile  occurs  at  slightly  above  100  kms  for  both  75  (Fig.  22)  and  45  Hz 
(Fig.  23). 

In  a narrow  angular  band  about  the  magnetic  equator  where  quasi 
transverse  propagation  applies  in  the  ionosphere  the  ambient  profiles 
behave  like  metallic  reflectors  and  the  re-radiation  theory  gives  very 
clearcut  pictures  of  the  reflective  process  in  terms  of  levels  and  half 
widths  (Figs.  39-41).  This  is  also  true  for  highly  depressed  ionospheres 
(Figs.  46-48). 
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APPENDIX  A 


Low  Altitude  Absorption 

The  low  altitude  absorption  occurs  where  the  electric  field  com- 
ponents of  the  rf  wave  satisfy  the  following  inequality 

lE3l>>  IE,  I,  IE2I  A.l 

where  the  indices  1,  2,  3 stand  for  x,  y and  z respectively.  Thus  the  low 
altitude  absorption  is  due  principally  to  the  vertical  electric  field 

component  so  that 

J Re  ( i E3  E3  M33  1 


Re 


.*  (1+M33>*  (I+M33)  33 


A.2 


M33)*  (I 

where  represents  the  vertical  electric  field  canponent  at  the  base  of 
the  ionosphere  (Ej^  to  very  good  approximation  may  be  taken  equal  to  the 
vertical  electric  field  at  the  ground).  The  substitution  of  Ejj/(l  + M33) 
for  E3  is  a cor  sequence  of  continuity  of  the  vertical  component  of  the 
displacement  vector  along  with  the  inequality  A.i.  We  will  first  limit 


the  discussion  to  the  case  for  which  u>  for  all  species.  Then 

N a J 

Mjj”'  -i  E Xj 
J=1 

where  J is  a species  index  including  both  electrons  and  ions  and 

X.  = ^ 

' "•'J 

Observe  that  X 4 is  real.  Substituting  A. 3 into  A.2  yields, 
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This  Dtaximizes  at  that  altitude,  z#,,  for  which 

E W"' 

j=i 

From  (A.5)  the  line  shape  in  the  neighborhood  of  Zq  may  be  written  as 


A.6 
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(normal iiation  1s  to  unity  at  2=2.) 

I 0 


ui 
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r N 1 

2 

L j=i  J 

where  the  prime  denotes  the  derivative  with  respect  to  2.  If  the  profile 
densities  and  collision  frequencies  are  assumed  to  vary  exponentially  in 


the  neighborhood  of  2^  with  scale  heights  that  are  invariant  over  the 


line  profile  then  A. 7 becomes 
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where  a and  y are  inverse  scale  heights  at  for  the  profiles  and  col 
lislon  frequencies  respectively.  Equation  A. 8 has  been  compared  with 
the  full  wave  low  altitude  absorption  line  shape  in  section  4. 

Equation  (5)  can  be  used,  along  with  the  preceding  equations  to 
estimate  tiie  integrated  effect  of  the  low  altitude  absorption  line, 
Since  decreases  quite  rapidly  above  the  level  and  since  it  is 
nearly  constant  between  z,^  and  the  ground  it  that 
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Using  the  Umn-zian  line  profiles  just  discussed^  the  integral  of,  J 
over  the ■■  line  becomes  - V'-  ^ >'v^ 
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if  A.7  is  used 
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Combining  factors,  Eq  (8)  yields  the  following  estimate  for  the  integrated 

low  altitude  absorption  (measured  in  terms  of  attenuation) 

^ 0.2flS|2  A n 

Re(S)/3zo  ‘^®^lOOOkm 

with  f the  frequency  in  Hz. 

Recall  that  the  preceding  equations  have  been  developed  subject  to 
the  condition  that  foif*  all  species.  This  condition  is  satisfied 

for  the  daytime  low  altitude  absorption  and  is  also  true  for  all  species, 
except  possibly  electrons,  for  the  low  altitude  nighttime  absorption.  If 
for  a given  species 

lnlcu^;j»Vj  A.U 

then  the  only  change  in  the  preceding  development  is  that  the  term  in 

p 

A. 3 corresponding  to  that  species  is  to  be  multiplied  by  n where  n is  the 
direction  cosine  between  the  geomagnetic  field  and  the  z axis.  .Although 
A. 14  is  only  marginally  satisfied,  this  factor  has  been  included  for 
electrons  in  the  45  and  75  Hz  I'esults  for  the  nighttime  profile  discussed 
in  section  4. 


APPENDIX  B 


High  Altitude  Absorption 

To  arrive  at  an  estimate  for  the  location  of  the  high  altitude 

absorption  and  a formula  for  its  line  shape,  it  will  be  assumed  that  it 

occurs  in  an  altitude  region  where  only  the  outgoing  whistler  mode  has 

survived.  The  index  of  refraction,  assuming  just  a single  species  of 

+ ions,  for  the  whistler  mode  is  (in  this  appendix  the  subscript  e is 

for  electrons  and  i for  ions): 
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Now  in  the  altitude  region  of  interest 
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and  to  a reasonable  approx iination 
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Thus  B.3  is  further  approximated  to  give 
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The  second  term  in  the  parenthesis  of  B.6  can  be  no  larger  than  0.5, 
thus  we  assume 
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and  finally  obtain 
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From  (B.6)  the  conductivity  may  be  identified  to  be 
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whereas  the  imaginary  part  of  N from  (B.8)  is 


111)  N = -r 
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Thus  in  the  high  altitude  region  the  absorption  per  unit  volume  varies  as 

, vCT;  wi  u>  ,»>i 
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(o)  + Slj)"  + vf 

The  problem  of  finding  the  height,  Zq,  where  the  high  altitude  absorption 
maximizes  is  tantamount  to  maximizing  the  function 
1 ^ \/u>Tr 

G(z)  = lii)l<*,2,;j/((w  + aj)2+y2)j ^ 7 Ptj 

*'  (tu +Ilj)  + i>j2 

Assuming  that  the  ion  density  and  collision  frequency  vary  exponentially 
in  the  neighborhood  of  z^  with  scale  heights  that  are  invariant  over  the 
line  profile  then  (0.10)  maximizes  at  that  height,  Zq,  for  which 


B.12 


L + CKj  - 7.  = 0 

^ (w  + J2j)^  + vf  (ca)  + + vf 

where  a and  y are  inverse  scale  heights  for  the  density  and  collision 


frequency  respectively. 

In  the  neighborhood  of  the  absorption  line  assumes  a Gaussian 
line  shape  with  half  width  determined  by  evaluating  the  second  derivative 
of  6 (i.e.,  B.ll)  at  Zq.  Carrying  through  this  procedure  yields  for  the 


high  altitude  line  shape 
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Equation  B.13  along  with  B.12  has  been  used  in  section  4 to  compare  with 
the  location  and  shape  of  the  full  wave  high  altitude  absorption.  The 
integrated  intensity  of  the  line,  unlike  the  low  altitude  absorption 
line,  cannot  be  so  easily  approximated  since  it  depends  upon  the  trans- 
mission through  the  0 and  E layers. 


APPENDIX  C 


Resonance  For  A Simplified  Model 

The  discussion  in  this  appendix  is  based  on  the  simplified  model 
shown  in  Fig.  (29).  It  will  be  assumed  that  the  altitude  range  - Z2» 
is  sufficiently  large  that  although  both  extraordinary  and  ordinary 
wave  exists  at  the  altitude  only  extraordinary  wave  exists  at  z-^- 
Denoting  the  extraordinary  wave  by  the  superscript  x and  the  ordinary 
wave  by  the  superscript  o,  the  total  rf  field  components  E^,  E^,  3C^,  3(’2 
at  the  ground  2=0,  in  terms  of  the  extraordinary  and  ordinary  field 
components,  Ej’®,  £2*°,  3^2*°*  ^1 

B|  = Hj  coskCzj  + B”  coskCzj  + slnkCZj  + sin  kC/j  c.i 

B.,  = VJ\  coska, +U“  coskCz, sin  kCz,  - JCj  sinkC/.,  c.2 

3f|  « 3f J coskQ, +3C‘J  coskCz, -iCB*  sinkCz, -iCB“  sin  kC/.,  c.3 

^2  *305  coskCz, +3d^  coskCV, ^ B*  sin  kC/.,  + ^ Bf  sin  kCz,  C.4 

With  allowance  for  co-ordinate  system  differences,  these  equations  can 
be  derived  following  Booker, 

In  the  following  development  it  is  assumed  that  the  extraordinary 
index  of  refraction,  nj^j,  and  the  ordinary  index  of  refraction  Ogj,  for 
regions  js2,  3 and  4 satisfy  the  condition 

In2jl,ln“jl>>(s2|  j»2.3,4 

Assuming,  then,  vertical  propagation  of  the  extraordinary  wave  between  Zj 
and  22  and  the  total  decay  of  the  ordinary  wave  at  23  the  field  com- 
ponents lij'*’  , li*'“  . 3f^*®  , aUj  way  all  be  expressed  in  terms 
of  two  independent  quantities  A and  B as  follows 
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E*  = i — (1-Q)  , E®  = -i-?- 

2 n^2  2 0^2 

C.6 

• 

3C*  = -iA(l  +Q)  , 3(0  =iB 

C.7 

» ; 

^2  = A(1  tQ)  , 3(0  =B 

C.8 

1 

where 

Q = R exp  (-2  ikn^2  ^ ^ 

C.9 

' 

r23  + r34  exp  (-2  ikn^j  (Z3  - Z2) ) 

Ra— 

1 tr23r34  exp  (-2  ikn^3  (Z3  - Z2) ) 

C.IO 

•I 

■'i 

"x3-'’x2 

r,,  0——— 

^ '\3  + “x2 

C.ll 

» 1 

»x4-«x3 

• 

■■i 

^34  ° ^ 

^ «x4  + '‘x3 

C.12 

'.i 

i 

Substituting  C.5  through  C.8  Into  C.l  through  C.4  yields 

t;  1“Q  ^ 1 * . f 1 . . 


fi“Q 

li,  “ — —coskC 
’ ih2 


coskCzj  fiC(l  4’Q)sinkCz 


1*4 


coskCZ]  + iCsiukCz,  B C.13 


^2  « I kC2|  - sin kCZ| j A +1“  ^cos kCzj  ♦ i sin kCzJ  B C . 14 

• x2  ^ 

a|~i(lfQ)uMkCZ|  + ^ (loQlsinkCzij  At |ii;o9kCZ| -••^~sinkCZ||  B 


3C  = 1(1 +Q)coskCz, +— ^(l-Q)sinkCzjl  A+  coskCzj+— !— sinkCzj  B 

2 [ Cn^2  J '‘”o2  J 

Assuming  a perfectly  conducting  ground,  the  conditions 

C.17 

must  be  fulfilled.  After  some  algebra,  the  conditions  C.17  can  be 
manipulated  into  the  following  form  for  the  modal  equation 

/2icoskCz,  sinkC/.  / I \\  / coskCV.|  / 1 X 

y - — r . — (c.  J-)) 

(2  i cos  k('z.  sin  kC/..  , i A 1 coskrZ|  , fX 

♦ (t’+  — j]  +sin  kC’zi  ( 2 i sin  kC/.i  + (C+  — j) 

''o2'‘x2  '‘x2  V ^7  \ ''o2  ' 

+ 1=0  C.18 


Further  reasonable  approximations  for  the  cases  of  interest  here  are 
kCz,  <<  I 

Using  these  conditions  it  is  a straightforward  exercise  to  show  that 
correct  to  order  (kz^  lnJ^2l)”^  (C.18)  becwiies 

Q [1  + -J — (i+r2)  Li«o 

I ^«k2  j 

or  equivalently  to  order  {kZj 
Qc,i>f.pj! — (i+c^)  1+1-0 

J 

Substituting  C.9  for  Q into  C.21,  the  latter  becomes 
oxp  I -2  iknjj2  (Z3  - 2| ) + ( 1 + C“^]+ KnR  j « -1 


Taking  C.22  yields  for  the  square  of  the  cosine  of  the 

eigenangle 


V 


k?.,nx2[  2kn^2(22-z,)-(2p+l)jr-  C.23 

where  p = 0,  1,  2 . . . 

If 

<<  1 and6<<  1 . C.24 

I Re (0x2)1 

2 

Eq.  C.23  yields  a sharp  resonance  behavior  for  C as  well  as  for  the 
attenuation, 

a=  -8.68kRo(S)X  10^  dB/IOOOkm,  C.25 


where  S is  the  sine  of  the  eigenangle.  If  in  addition  to  the  conditions 
C.19. 

C.26 


the  resonance  condition  is 

I / \ 

4 (Re(..,2)j 


Xff  free  space  wuvcicagtii 


C.27 


Since  AQ/Re(nj^2)  the  local  wavelength  in  region  2,  a primary  condition 
for  resonance  in  this  instance  is  that  region  2 serve  as  a quarter  wave- 
plate.  The  first  condition  of  C.24  is  always  true  for  cases  of  interest 
here  so  that  the  second  condition  of  C.24  becoaies  the  crucial  one.  For 
the  cases  of  interest  here*  namely 

Injj3l>>  In^jl,  Inx4l  C.28 

■94^ 


the  second  condition  of  C.24  will  be  best  satisfied  when  the  layer, 
region  3,  also  has  a thickness  close  to  quarter  integer  multiples  of  the 
local  wavelength  in  region  3,  Also,  the  condition  most  detrimental  to  a 
sharp  resonance  is  the  condition  that  the  layer,  region  3,  has  a thick- 
ness close  to  half  integer  multiples  of  the  local  wavelength.  To  under- 
stand better  the  origin  of  these  conditions,  consider  C.IO  subject  to 
the  qualification  that  regions 2 and  4 have  identical  charge  densities 
and  collision  frequencies.  Then  subject  to  the  first  condition  of  C.24 
and  conditions  C.28  it  follows  for 

-2  ikrtjjj  (Z3  - Z2)  = (-2  pir -e)  i - 2k  1 Iiu  (njjj)l(z3  - ^2^  • p=  1.2,3...  C.29 
and  for 

e,  2k  1 Im  (nj^3 ) 1 (Z3  - Z2 )< < I 
that  C.IO  becomes  approximately 

2 kl  Im  “ *->)  + ie 

g Q,  — :: ^ 

2 kllm 

For  conditions  of  primary  Interest  here 

2 ki  1111(11^3) Uz) 

If  in  addition 

c<<4ln^2\3l  C.33 

it  follows  that  R<*:1  so  that  the  6 term  in  C.23  will  tend  to  smear  the 
resonance. 

If  instead  of  C.29  the  condition 

“2  lkii^3  (Z3  - Z2 ) *“  -( ( 2p  + 1 ) w + e)  i - 2 k I Im  (n^s ) UZ3  - z, ^ • p “ 0, 1 ,2  . . 
is  imposed  and  conditions  C.30  satisfied  it  follows  that  R°1  to  first 

order  in  c,  kn|„(»\3)  1(23*22)  and  nj^g/n^^j  . Thus,  subject  to  these 
conditions,  the  resonance  peak  should  be  accentuated  when  C.34  is 
satisfied.  It  is  interesting  to  note  that  contrary  to  the  preceding 

95 


C.30 

C.31 

C.32 


findings,  the  condition  C.29  rather  than  C.34  was  previously  suggested 

25 

as  being  the  condition  for  strong  resonance.  In  section  4 results  of 
calculations  based  on  C.23  will  be  given  which  illustrate  the  foregoing 


behavior. 


For  the  more  general  case  of  two  reflecting  layers  separated  by 
an  inhomogeneous  medium,  we  would  expect  C.25  to  be  replaced  by 


Re  (iijj2)dz  = (2  p + I )X^,/'4  ; p = 0,  I,  2, 3,  etc. 


To  examine  the  polarization  of  the  resonant  wave,  it  is  first 
observed  that  C.20  gives  approximately 

0 = -l +i  (I +C’““)/kZ|n^2  C. 

Inserting  this  in  C.14  and  remembering  that  £2  = 0 gives 
A = - ikz,  C- 1«^2  B / S“  C. 

Using  C.37  in  Eqs.  C.15  and  C.16  and  retaining  just  the  dominant  tenns 
give 

3f|/3f2“ -ic:-  C. 

2 

Since  |C  1 is  on  the  order  of  unity  for  the  resonant  wave  for 
reasonable  values  of  the  paranieters,  Eq.  C.3B  shows  that  the  polariza- 
tion close  to  resonance,  unlike  the  normal  situation  at  ELF,  will  con- 
tain comparable  anounts  of  both  vertical  and  horizontal  polarization. 
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